Proposal to determine the Fermi-surface topology of a doped iron-based 
superconductor using bulk-sensitive Fourier-transform Compton scattering 
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We have carried out first-principles calculations of the Compton scattering spectra to demonstrate 
that the filling of the hole Fermi surface in LaOi-^F^FeAs produces a distinct signature in the 
Fourier transformed Compton spectrum when the momentum transfer vector lies along the [100] 
direction. We thus show how the critical concentration x c , where hole Fermi surface pieces are filled 
up and the superconductivity mediated by antiferromagnetic spin fluctuations is expected to be 
suppressed, can be obtained in a bulk-sensitive manner. 
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PACS numbers: 71.18.+y, 71.20.-b, 74.25.Jb 74.70.Dd 

The Fermi surface (FS) topology is a key ingredient for 
high temperature superconductivity in iron based lay- 
ered pnictides. The so-called s± modeli"— predicts su- 
perconducting gaps of one sign on the FS cylindrical hole 
sheets near T(0, 0) and of another sign on cylindrical elec- 
tron sheets at M(7r,7r). Doping x is needed to move the 
system away from the magnetic instabilities due to FS 
nesting^"—. In the superconducting material, spin fluc- 
tuations (related to residual FS nesting) may provide a 
glue for the Cooper pair si. However, a complete filling of 
the hole FS at a certain electron doping will eventually 
lead to the suppression of the spin fluctuation glue. The 
exact value x c of this critical doping might be affected 
by subtle correlations effects^. 

Experimental information regarding the FS topol- 
ogy comes mostly from angle resolved photoemission 
spectroscopy (ARPES) which is a surface sensitive 
probe . Moreover, since the doping level in the 
bulk could be different from that at the surface 2 -, the FS 
signal should be checked with bulk probes. Since the FS 
information from quantum oscillation studie s 13 ! 14 can be 
distorted because of the required high magnetic fields, 
we suggest determination of the FS topology via Comp- 
ton scattering measurements ] 15 i 16 In this letter, we show 
that a one dimensional Fourier transform of the Comp- 
ton profile along [100] presents a large signal when the 
hole Fermi surface vanishes, providing a bulk sensitive 
method for determining the critical doping x c for high 
temperature superconductivity in LaOi-^F^FeAs. 

Recent advances in synchrotron light sources and 
detector technology have renewed interest in high- 
resolution Compton scattering as a bulk probe of fermi- 
ology related issues, see e.g., RefsJ^"— . In a Comp- 
ton scattering experiment, one measures a directional 
Compton profile (CP), J(p z ), which is related to the 
twice integrated ground-state electron momentum den- 



sity p(p x ,Py,Pz) by 



J(Pz)=J J p(Px,Py,Pz)dp x dp y 



(1) 



for high momentum and energy transfer 21 . The explo- 
ration of FS topology with the aid of Compton scatter- 
ing is complicated by the double integral in Eq. [TJ As 
a result, FS breaks in pijPxi Py,p z ) do not usually induce 
rapid variations in J(p z )- A possible approach to deal 
with this problem is to measure CPs along many different 
directions and use state-of-the-art reconstruction meth- 
ods based on the autocorrelation function B(x,y,z) to 
obtain p(p x ,Py,Pz) 22 • A much simpler and more robust 
method will be proposed in this letter. 

The calculations presented here were performed within 
the local density approximation (LDA) using an all- 
electron fully charge self-consistent semi-relativistic 
Korringa-Kohn-Rostoker (KKR) method 23 . The com- 
pound LaOi-^F^FeAs has a simple tetragonal structure 
(space-group P4/nmm). We have used the experimen- 
tal lattice parameters 24 of LaO0.87F0.13Fe As in which no 
spin- density- wave order was observed in neutron scatter- 
ing experiments. Self-consistency was obtained for x = 
and the effects of doping x were treated within a rigid 
band model by shifting the Fermi energy to accommo- 
date the proper number of electrons 2 ^"—. The conver- 
gence of the crystal potential was approximately 10 -4 
Ry. The electron momentum density (EMD) p{p x ->Py>Pz) 
was computed on a fine mesh of 40.4 x 10 6 within a sphere 
of radius 12 a.u. in momentum space. To simulate the 
effect of resolution in high resolution Compton scattering 
experiments, the EMD has been been convoluted with a 
Gaussian characterized by a FWHM of 0.17 a.u. 

Since the LaOi-^F^FeAs electronic structure has a 
two-dimensional character, we shall focus on the calcu- 
lated (001) 2D-projection of the momentum density given 
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FIG. 1: (Color online) Calculated LCW distribution of in the 
paramagnetic tetragonal Brillouin zone of LaOFeAs. 



by 



p 2D (Px,Py) = J p(p x ,Py,p z )dp z . (2) 



In experiments, one adopts the so-called direct Fourier- 
transform method 22 to reconstruct p 2D (p x ,Py) from sev- 
eral directional CPs J(p z ) measured in the (001) plane. 
This method uses the autocorrelation function B(x, y, z) 
which is straightforwardly denned as the Fourier trans- 
formation of the momentum density 



B(x,y,z) = J J J dp x dp y dp z (3) 

x p(Px,Py,Pz)e*p[i(PxX+P y y+Pzz)] • 

Since p(p x ,Py,Pz) can be expressed as a sum over the 
momentum density of the natural orbitals ipj(x,y, z)2& 
by using the convolution theorem, it can be shown that 
B(x,y,z) is the autocorrelation of the natural orbitals 



B(x,y,z) 



?•*/// 



dudvdw 



(4) 



x ijjj(x + u,y + z + w)ip* (u,v,w) 



where nj is the occupation number of the natural orbital 
ijjj(x,y,z). In the experiments B(x,y,z) is obtained di- 
rectly along a given direction by taking the ID-Fourier 
transform of the CP along that direction. Then, once a 
set of £>'s has been calculated, a fine mesh is set up in 
real space and B(x,y,z) is obtained at every mesh point 
by interpolation. Finally, if desired, an inverse Fourier 
transform of B(x, y, 0) yields the distribution p 2D (p x ,Py)- 
Our simulations reveal that the breaks in p 2D (p x ,p y ) 
caused by FS crossings are scattered throughout momen- 
tum space with small weights given by matrix elements 
involving mostly the Fe d orbitals. Therefore, FS features 
are not easily detected directly in the p 2D (p x ,p y ) distri- 
bution. However, as shown in Fig.[T]the Lock- Crisp- West 
(LCW) folding 31 can enhance FS breaks by coherently 
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FIG. 2: (Color online) Ap\%l(p x ,p y ) for x 2 = 0.15 and x\ = 
0.10. The yellow labels indicate the high symmetry positions 
r and M in momentum space. 
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FIG. 3: (Color online) AB\%*(x,y) for x 2 = 0.15 and xi = 
0.10. The yellow labels Fe indicate the atoms in the iron 
sublattice; green labels As indicate the (x, ^-projection of 
the atoms in the arsenic sublattice. Atomic assignments are 
based on assumption Fe is at origin. 



superposing the umklapp terms according to 

n(k x ,k y )= P 2d ( k x + G x ,k y + G y ) , (5) 

where n(k Xl k y ) gives the number of occupied states at 
the point (k x ,k y ) in the first Brillouin zone by summing 
over all projected reciprocal lattice vectors (G x , G y ). The 
maximum of n(k Xl k y ) at M(7r,7r) is associated with the 
electron pockets while the minimum at V (0,0) is related 
to the hole pockets. Since the LCW folding can also en- 
large artificial errors from the experimental data, below 
we will introduce a more robust means of extracting in- 
formation about the evolution of the FS topology with 
doping. 
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FIG. 4: (Color online) AB\%*(x,y). for x 2 = 0.15 and xi = 
0.10, after artificially removing all the Bloch states around M. 
Atomic assignments are based on assumption Fe is at origin. 

We can get more precise information on wave function 
symmetry near the Fermi surface (FS) by taking differ- 
ence maps between two nearby dopings, 

&p\Zl(Px,Py) = p 2D (Px,P y )\x 2 - p 2D (p x ,p y )\ Xl , (6) 

where X2 and x\ are two different doping levels such that 
X2 > x\. The subtraction in Eq. [6] acts as a projector on 
the Fermi level subspace with the advantage of eliminat- 
ing the large isotropic contribution of the core and some 
irrelevant valence electrons. The difference for x^ = 0.15 
and x\ = 0.10 shown in Fig. [2] displays interesting FS 
effects strongly modulated by Fe d wave function effects. 
The corresponding Fourier transform AB\%1 (x, y), Fig.[3j 
separates the different length scales in real space, which 
contribute to the oscillations in Ap\%*(p x ,p y ). Thus, 
the peaks in the autocorrelation function AB indicate 
characteristic distances over which wave functions at the 
Fermi level are coherent. The peaks in Fig. [3] mostly 
stem from the Fe d orbitals since these largely dominate 
at the Fermi level^. In fact, from Fig. [3] one can see that 
the main peaks correlate very well with the iron sub- 
lattice. However, note that there are weaker features, 
marked "As", which correlate with the positions of the 
As atoms. Since the bands near the Fermi level mostly 
consist of Fe d xz , d yz and d x i_ y i orbitals 32 , Fig. [3] re veals 
these characters in real space. To facilitate comparisons 
we normalize AB to unity at the origin. 

We employ filtering techniques to enhance the sensi- 
tivity of the A£?-maps to particular FS cylinders. Thus, 
in Fig. [4] we artificially remove the FS around M by ap- 
plying a filter cutting out the Bloch states near M. As a 
result some Fe peaks essentially disappear, revealing the 
wavefunction characters of the T cylinders. In the same 
way, Fig. [5] shows the corresponding maps for the elec- 
tron cylinders, generated by filtering out the cylinders at 
T . This filtering procedure can be tested by calculat- 
ing AB\**(x,y) at higher doping for X2 and xi, where 
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FIG. 5: (Color online) AB\Z*(x,y). for x 2 = 0.15 and xi = 
0.10, after artificially removing all the Bloch states around T. 
Atomic assignments are based on assumption Fe is at origin. 



the number of added electrons is sufficient to remove the 
hole cylinders at T without any filtering needed* ^ 11 ! 33 
The result is very similar to Fig. [5j confirming that the 
present filter is an efficient way of sorting out contribu- 
tions from different FS cylinders 34 . By comparing Fig. [4] 
and Fig. [5j it is clear that the hole FSs and the electron 
FSs give strikingly different contributions to AB\%1 (x, y). 

In Fig. [4] [hole FSs], there is a peak at every Fe site, 
with a corresponding weaker network of peaks at As sites. 
In contrast, for the electron contribution in Fig. [5j there 
are dips at some Fe sites and the As sites signatures are 
not visible. Some of these features can be readily un- 
derstood with the phase factor exp(zk • r) of the Bloch 
wave function. At the band bottom, T, k = and all 
atoms are in phase (bonding). At the top of the band, 
(7r,7r), the phase factor is (— l) m + n - i.e., perfectly an- 
tibonding. Since the FSs are small cylinders near T and 
(-7T, 7r) respectively, their A£>-maps are dominated by this 
interference term. 

Interestingly, in Fig. [6j we show 35 the contrasting be- 
havior of the autocorrelation function AB near the sec- 
ond neighbor Fe, at a distance about 7.6 a.u. from the 
origin along the [100]-direction. The minimum is pro- 
duced by anti-bonding states belonging to the M FSs 
while the maximum is the result of bonding states be- 
longing to the T FSs. Hence with doping the feature 
should evolve from the red curve when both FS are oc- 
cupied to the green curve when the T FS disappears. 

Since this last analysis will need only one FT of the 
experimental Compton profile along [100] it clearly pro- 
vides particularly robust FS information in real space. 
In contrast, the p 2D reconstruction requires many more 
numerical manipulations 2 - and therefore it can be less 
reliable. Our AB methodology thus is a more reliable 
way to detect when the FS signal at T vanishes. 

In conclusion, our study predicts that high-resolution 
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FIG. 6: (Color online) Cut of AB along the [100] direction. 



Compton scattering spectra will yield signatures for the 
FS topology of Fe-based superconductors. In particular, 
the AB-map projected along the [100] direction displays 
a remarkable signature of the FS evolution with dop- 



ing. Thus our method gives a robust way to establish 
the topology rather than the precise shape of the iron 
pnictides FS. These results indicate that Compton scat- 
tering can provide a powerful new spectroscopic window 
for investigating FSs compatible with the s± model for 
the description of the superconducting order parameter 
in the Fe-based superconductors. 
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